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Abstract 

The quark Dyson-Schwinger equation shows that there are distinct differences between 
light and heavy quarks. The dynamical mass function of the light quarks is characterised 
by a sharp increase below 1 GeV, whereas the mass function of the heavy quarks is 
approximately constant in this infrared region. As a consequence, the heavy-meson masses 
increase linearly with the current-quark masses, whereas the light pseudoscalar meson 
^ masses are proportional to the square root of the current-quark masses. 

Quark propagator 

The Dyson-Schwinger equations [DSEs] provide an excellent tool to study nonpertur- 
bative aspects of quark propagators and bound states 1 . Various properties of the light 
mesons can be described very well in this framework 2 but this method is not limited to 
^ ■ the light quark sector. Nonperturbative effects, such as chiral symmetry breaking and 
confinement, are reflected in the dressed quark propagator S(p)^ 1 = z 7 ■ pA(p 2 ) + B(p 2 ). 
In Euclidean space, the DSE for the quark propagator is 

>] Sip)' 1 = Z 2 ij ■ p + Z.mifi) + Z x J ^g 2 D, u (p-q)^S(q)T a u (q,p), (1) 



On 



(N 

\o 

where D fll/ (k) is the dressed gluon propagator, T®(q;p) is the dressed quark-gluon vertex 
and m(fi) is the renormalised current-quark mass. The quark-gluon-vertex, quark wave- 
function, and quark mass renormalisation constants, Z\, Z 2 , and Z 4 respectively, depend 
on the renormalisation point //, as does the current-quark mass. 

The set of DSEs form an infinite hierarchy of integral equations, and for actual calculations 
we need to truncate this set of equations. Following Ref. 1, we use a truncation based on 
^ \ the so-called ladder approximation, namely 

Zv K(q,p)g 2 D^(p - g)f 7 , -> £ 7 „0((p - q) 2 ) D^(p - g)f 7 , . (2) 
• 1— 1 . 

^ ■ Qualitatively reliable studies of the DSE for the dressed gluon propagator indicate that 
the effective interaction is significantly enhanced in the infrared such that on this do- 
main it is well represented by an integrable singularity 3 . Combining this observation 
with asymptotic freedom in the ultraviolet, motivates an Ansatz 2 for Q((p — q) 2 ) with a 
8— function to represent the integrable infrared singularity, plus a tail term, regular in the 
infrared region, to incorporate asymptotic freedom. The renormalisation constants are 
fixed by applying the renormalisation condition A(fi 2 ) = 1 and B(fi 2 ) = m(/i). 
Now that we have fully specified the truncation of the quark DSE, we can solve numer- 
ically the resulting coupled integral equations for A(p 2 ) and B(p 2 ). In Fig. 1, we show 
the dynamical quark mass function M(p 2 ) = B(p 2 )/A(p 2 ) obtained for different quark 
flavours. With a nonzero current-quark mass, the dynamical mass function behaves as 

M(p 2 ) * r - 1x7m , P 2 > Aqcd , (3) 
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where 7m = 33 _ 2 ^ f is the anomalous mass dimension, and m the renormalisation-point- 
independent current-quark mass. Dynamical chiral symmetry breaking is manifest when, 
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Fig. 1. The dynamical (ffiskfib function of the quark propagator in the qhj©al/)limit (circles), and 
for up/down (solid), strange (short dash), charm (dash-dot), and bottom (long dash) quarks. 



for m = 0, one obtains M(p 2 ) ^ in solving Eq. (Jj), which is impossible at any finite 
order in perturbation theory. Analysing our chiral limit solution we find 
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P 2 > Aqcd 



(4) 



with the renormalisation-point-independent chiral condensate — (qq)° = (0.227 GeV) 3 . 
The large-p 2 momentum dependence of the quark mass function, both Eq. (H) and Eq. (§]), 
is characteristic of the QCD renormalisation group at one-loop 4 . 

In the infrared region, we see that dynamical chiral symmetry breaking implies more than 
just (qq)fj, 7^ 0. For small momenta, p < 1 GeV, there is a qualitative and distinct dif- 
ference between light and heavy quarks. The light-quark mass function is very similar 
to the chiral limit mass function, and is characterised by a significant infrared enhance- 
ment, a direct result of that in the effective coupling, Q(k 2 ). Introducing the Euclidean 
constituent- quark mass, M E , as the solution of p 2 = M 2 (p 2 ), the ratio: M E /m(/i), is a sin- 
gle, indicative and quantitative measure of the nonperturbative effects of gluon- dressing on 
the quark propagator. For current-quark mass values of m u /<j = 3.7 MeV, m s = 82MeV, 
m c = 0.59 GeV and mj = 2.0 GeV at the renormalisation point fi = 20 GeV, we find 
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which clearly indicates the magnitude of this effect for light quarks. It provides a qualita- 
tive understanding of the nature of the "constituent-quark" mass. The magnitude of this 
ratio is a signal of the effect of the infrared enhancement of the quark-quark interaction 
on quark propagation characteristics. 

For heavy quarks, the ratio M E /m(fi) takes a value of 0(1) because the current-quark 
mass is much larger than the mass-scale characterising the infrared enhancement in 
the effective coupling, Aqcd- This means that in the spacelike region the momentum- 
dependence of the heavy-quark mass function is dominated by perturbative effects and 
suggests, in phenomenological applications, that one can replace the heavy-quark mass 
function by a constant mass 5 . 
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Note that the both heavy and light quarks, and also the gluons, are explicitly confined; 
the essence of confinement is that the propagator of a confined particle does not have 
a Lehmann representation 6 . Therefore, quark confinement entails that there is no pole- 
mass, which would be the solution of p 2 + M 2 (p 2 ) = 0. 

Meson bound states 

The renormalised Bethe-Salpeter Equation [BSE] for meson bound states is 

[lVp;P)L = J-0^[S(q + )r H ( q] P)S(q^] sr K^( q ,p-,P), (5) 

where: H specifies the flavour structure of the meson; q + = q + r\ P P, g_ = q — (1 — 7] P ) P, 
with P the total momentum of the bound state and rjp the momentum partitioning; and 
r,. . . ,u represent colour-, Dirac- and flavour-matrix indices. Solutions of Eq. (|^) exist only 
for particular, separated, timelike values of P 2 \ and the eigenvector associated with each 
eigenvalue, the Bethe-Salpeter amplitude [BSA]: Th{p',P), is the one-particle-irreducible, 
fully-amputated quark-meson vertex. Once we know this BSA, we can use it to calculate 
meson observables such as the decay constant and various form factors and other physical 
observables. For pseudoscalar mesons, the solution of Eq. (H) has the general form 2 

T H (q; P) = lh [iE H { q] P) + 7 • P F H (q; P) + 1 ■ qG H (q; P) + a^P u H H (q; P)\ . (6) 

In order to calculate the BSA, we must have an explicit form for the kernel K^(q,p; P) 
in Eq. fl5|). Requiring that the axial- vector Ward indentityQ be preserved constrains this 
kernel. The ladder approximation fulfills this constraint 7 

r H (p; P) = -\ J g((p - qf) Dj?(p - qh,S(q + )T H (q; P)S{q^) lu . (7) 

Numerical solutions of Eq. (0) for increasing values of the current-quark mass show that 
the meson mass grows with the quark mass. In Fig. 2 we present the meson mass m# as 
function of the renormalisation-point-independent current-quark mass rh, for both (uq) 
(ii ± , K, D, B) and (fictitious) (qq) (tt°) pseudoscalar mesons. In both cases, the meson 
mass mn is for light quarks, proportional to the square root of the current-quark mass, 
whereas for heavy quarks, the meson mass grows linearly with the current-quark mass: 
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rh q , m q < 0.15 GeV ~ m s , 
m q , m q > 0.3 GeV ~ 2m s 



The pseudoscalar decay constant increases with quark masses, at least up to quark masses 
of about 0.5 GeV ~ 3m s . The reversal of this trend, that is dictated by the heavy quark 
limit ju\/ m H — constant, 5 is reached at larger current-quark masses. 
Although the BSA does depend on the momentum partitioning rjp in Eq. (19), the physical 
observables, such as the meson mass and decay constant, do not, provided one takes into 
account the nonleading Dirac amplitudes Fh and Gh in Eq. ([]). These amplitudes are 
also essential in order to satisfy the axial-vector Ward identity, and are quantitatively 
important for physical observables 2 . 



*In the chiral limit, this Ward identity relates the meson BSA to the quark propagator: the pions are 
the Goldstone bosons associated with dynamical chiral symmetry breaking. 
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Fig. 2. The pseudoscalar meson 
mass for (uq) (plusses) and equal 
quark (qq) (diamonds) bound 
states, as a function of the current- 
quark mass m q . 
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The ultraviolet asymptotic behaviour of all amplitudes E(p;P), F(p;P), p 2 G(p;P) and 
p 2 H(p;P), is fully determined by the chiral-limit quark mass function: all behave like 
Eq. (Q), both for light and heavy pseudoscalar mesons. A difference between the BSA of 
heavy mesons and light mesons is that heavy mesons are broader in momentum space: 
they approach the asymptotic behavior at larger momenta, and are thus narrower in 
coordinate space than light mesons. 

Conclusion 

The dynamical mass function is a useful measure of nonperturbative effects. In the in- 
frared region, the mass functions of the light flavours, in particular the u- and d-quarks, 
are very similar to the chiral limit mass function M Q (p 2 ). They are all characterised by a 
sharp increase below 1 GeV: the constituent quark mass is much larger than the current- 
quark mass. The mass function of the heavy quarks does not show such behaviour: 
nonperturbative effects are much less important for the heavy quarks. For the 6-quarks, 
and to a lesser extent the c-quarks, the mass function can be approximated well by a con- 
stant: the constituent mass is approximately the same as the current-quark mass. One of 
the consequences of this difference is reflected in the mass spectrum of the pseudoscalar 
mesons: the heavy meson masses increase linearly with the current-quark masses, whereas 
the light meson masses are proportional to the square root of the current-quark masses. 
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